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We present a field-theoretical treatment of the critical behavior of three-dimensional weakly diluted 
quenched Ising model. To this end we analyse in the replica limit n —> the 5-loop renormalization 
group functions of the 4 -theory with 0(n)-symmetric and cubic interactions (H. Klemert and V. 
Schulte-Frohlinde, Phys. Lett. B 342, 284 (1995)). The minimal subtraction scheme allows one to 
develop either the y^e-expansion series or to proceed within the 3d approach, performing expansions 
in terms of renormalized couplings. Doing so, we compare both perturbation approaches and discuss 
their convergence and possible Borel summability. To study the crossover effect we calculate the 
\Q . effective critical exponents. We report resummed numerical values for the effective and asymptotic 

i-H ' critical exponents. The results obtained within the 3d approach agree pretty well with recent Monte 

Carlo simulations, y^-expansion does not allow reliable estimates for d = 3. 
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I. INTRODUCTION 



On 
O 
0\ 

g\ 

Influence of a weak quenched disorder on magnetic second order phase transition remains the subject of much interest. 
One of the central problems here is the answer to the questions: (i) do the critical exponents (of a homogeneous magnet) 
change under dilution (by a nonmagnetic component)? and, if yes, (ii) how do they change? 

I 

a 

[ causes changes in the critical exponents. This statement is known as Harris criterion. Later, for a large class of 
• • ■ d-dimensional disordered systems it was proven that the correlation length critical exponent v must satisfy the bound 
v > 2/d 0|. Both statements focus attention towards studies of the d = 3 Ising model, where the typical numerical 
i values of the above exponents, together with the magnetic susceptibility and the order parameter critical exponents 
J_! 1 in the pure case read: a = 0.109 ± 0.004, v = 0.6304 ± 0.0013, 7 = 1.2396 ± 0.0013, = 0.3258 ± 0.0014 lCl 

a ; 

The reply on the second question concerning the numerical values of the critical exponents of d = 3 weakly diluted 
quenched Ising model (random Ising model - RIM) is more complicated. Below we briefly review some experimental 
P-pT| theoretical |22}-|37|| and numerical [ 38 - 5lJ results for the critical exponents of d — 3 magnets described by the 



Replying the first question it was argued Q almost 25 years ago that if the heat capacity critical exponent a of the 
pure (undiluted) system is positive, i.e. the heat capacity diverges at the critical point, then a quenched disorder 



RIM. Although the statements of Refs. [l]^2 should hold in principle for arbitrary weak disorder, the new critical 
exponents appear only in a region of temperatures controlled by the concentration of the nonmagnetic component. 
Such a region not always may be reached in practice, where only the effective exponents are observed |5Sq| . Theoretical 
results are due to the application of the renormalization group approach. These are numerous for the asymptotic 
values of critical exponents but quite seldom for the effective ones J37) . 



This paper has been motivated by recent Monte Carlo calculations of the asymptotic |5lj] and effective |46j,f49|,|5C|] 
critical exponents for the three dimensional Ising systems with quenched disorder. Here, we deal with a renormalization 
group study of both asymptotic and effective critical behavior. Perturbation theory series, which are known in the 
5-loop approximation for the 0(n) model with cubic anisotropy [ p3| enable us to perform our studies in the replica 
limit n — > 0. The minimal subtraction scheme ]54[ ] allows to develop either the y^-expansion or to proceed in the 
3d approach, performing expansions in terms of renormalized couplings. Doing so, we compare both perturbation 
approaches and discuss their convergence and possible Borel summability. 



1 



The paper is arran ged as follows. In section || we give a brief account of some results on critical properties of systems 
of interest, section II] describes the model and the renormalization procedure. In section IV we analyze the series 
obtained, perform their resummation, and give results for the asymptotic values of critical exponents. Effective critical 
behavior is discussed in section [v| Section |v| concludes our study. Details of the resummation procedures exploited 
in our study are given in the Appendix A. 



II. REVIEW 



Experiments. Crystalline mixtures of two compounds provide a typical experimental realizations of a RIM (see Table 
Q). The first compound is an "Ising-like" anisotropic uniaxial antiferromagnet with dominating short-range interaction 
(e.g. FeF2, MnFz), the second one [ZnF-i) is nonmagnetic. Mixed crystals (Fe x Zni- x F2, Mn x Zn\- x F-i) can be 
grown with high crystalline quality and very small concentration gradients providing an excellent realization of random 
substitutional disorder of magnetic ions {Fe +2 , Mn +2 ) by nonmagnetic ones (Zn +2 ). The experimental evidence of 
the new critical behavior at a weak quenched dilution was the nuclear magnetic resonance measurement of the 
magnetization in Mjiq. 864-^^0. 136-^2 [[§• The value of the magnetisation exponent was found to differ strongly 
from that in undiluted sample (see Table || for details). In a few years this result was corroborated by nuclear 
scattering measurements of magnetic susceptibility and correlation length critical exponents in Fe x Zn\^ x Fi J6|,[l^] 
and Mn x Zn\- x F2 R| at different dilutions 1 — x. The linear birefringence measurements brought about the cusp-like 
behavior of a specific heat at the transition point with a = —0.09 ±0.03 In particular this proved that within the 
error bars the hyperscaling relation dv + a = 2 is satisfied. 

Crossover phenomena in diluted systems are governed in addition to the temperature by the concentration of the 
magnetic component. The experimentally obtained exponents are often reported to be effective ones (i.e. temperature- 
and dilution-dependent). However, already in the first experiments on the critical behavior of quenched diluted Ising- 
like antifcrromagnets it appeared possible to reach the asymptotic region. Thus, studying the critical regime in Ref. jq| 
the authors found neither a region in reduced temperature r where one finds "pure" Ising exponents nor any evidence 
of crossover from pure to random exponents. This was explained by the crossover either taking place outside the 
critical region or being too slow. In Ref. |Tl| l the crossover from pure to diluted critical behavior was studied and a 
magnetization exponent (3 which does not change under dilution for x < 0.05 was found. The crossover occurs within 
a very narrow range of r at relatively large values of r. In Ref. excellent agreement of the measured exponents 
7 and v with the theoretical values of the RIM was obtained for the temperature range 4 • 10~ 4 < r < 10 -1 . The 
value of the critical exponent /3 was also the subject of a crossover analysis in Refs. Ref- |§) concludes that the 

experimental errors are too large in order to distinguish between the pure Ising model and the RIM critical behavior. 
In Ref. (h]] no crossover was found after the correction-to-scaling has been taken into account, and the RIM critical 
behavior was found in the whole temperature range. 

It is known that the diluted Ising magnet in a uniform magnetic field H along the uniaxial direction exhibits static 
critical behavior of the random-field Ising model |55|, [3(|. The random-field Ising model is the subject of intensive 
recent experimental studies (see |57j for a review). Such experiments give additional information about the critical 
behavior of the RIM when performed for H = @-|||. 

Renormalization group results. The change of the Ising-type magnets critical exponents under dilution found its 
theoretical confirmation in the renormalization group (RG) calculations. There, the change of the universal properties 
is interpreted as a crossover from the "pure" fixed point present in the undiluted model towards the "random" one 
which characterizes the new critical behavior of the RIM. The corresponding RG equations are de gen erated on the one 
loop level, which leads instead of the familiar expansion in e = 4 — d to the expansion in ^fe [p2|-pit . This expansion 
was recently extended from the 3-loop |^,|26) to the 5-loop level |27],|8| for the critical exponents. The y / e : -expansion 
series for the amplitude ratios |5£| is available with the 3-loop accuracy as well (5l|. Although the -^/e-expansion 
series are known with pretty high accuracy they seem to be of no use for the d — 3 case |2{|. Alternatively, the 
scaling field RG approach lead (via the two fixed points scenario) to the d — 3 critical exponents values v = 0.70, 
7 = 1.38 |30| definitively different from their "pure" counterparts. 

RG equations of the massive d = 3 field theory [[30| appeared to be the most fruitful tool takling the problem. The 
resulting expansions, considered as the asymptotic ones (see, however, Refs. |2£,61, 32| and section IV of this paper) 



were resummed succesively increasing the order of the perturbation theory series from the 2-loop (2LA) |3lf |32[| 
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through 3-loop (3LA) |53|,[54j to the 4-loop (4LA) approximations |3f],|36|]. Depending on the resummation procedure 
applied, the results for the correlation length and the magnetic susceptibility critical exponents {is, 7} read: 2LA: 
{0.67813, 1.33551} @ {0.679, 1.337} f|; 3LA: {0.670, 1.325} §§; {0.671, 1.328} @; 4LA: {0.6702, 1.3262} [||; 
{0.6680, 1.318}, {0.6714, 1.321} |3(|. Critical amplitudes universal ratios at d = 3 were obtained in 3-loop |]3| and 
5-loop approximations. 

It is worth to mention here that the renormalization at fixed d = 3 dimensions does not necessarily mean an implemen- 
tation of the massive scheme |30| in which most of the theoretical results |3l]-|36| were obtained. Indeed, the minimal 
subtraction scheme |54| is also suited for the d = 3 renormalization and can be applied without e-expansion ||S5| . 
Both schemes serve as complementary tools to obtain information about the critical behavior of the pure Ising model 
(or, more generally, 0(n)-symmetric </> 4 theory). However in the case of the RIM there exists only one study relying 
on the 3d RG approach in minimal subtraction scheme. In particular the resummed values of the critical exponents 
obtained there in a 3-loop approximation read: v = 0.666, 7 = 1.313 J3?j. 

Griffiths singularities and replica symmetry breaking. Although the RG approach does not allow to determine the 
region of concentration where the scenario with two fixed points, explained above, is applicable, it is generally 
assumed that it holds at least for weak dilutions. The idea of RG presumes that one studies the influence of (thermal) 
fluctuations around the spatially homogeneous unique ground state. This holds for the pure system but does not hold 
for the diluted one. Here, in a disorder dominated region one finds a macroscopic number of spatially inhomogeneous 
ground states, corresponding to the local minimum solutions of the saddle-point equation for an effective Lagrangian 
|p6|,|6^|. Physically the inhomogenous ground state corresponds to the so-called Griffiths phase |38| caused by the 
existence of ferromagnetically ordered "islands" in the region of temperatures between the critical temperatures of 
pure and diluted systems. This is described by a replica-symmetry breaking Lagrangian [36| leading to a behavior 
at the critical point (6^, [70| different from the pure one. The "traditional" RG results [^2j-37,[58|,|5!^,[6^,Q are valid 
for the replica-symmetric Lagrangians |7l| . Since the usual RG approach is based on an integration over the disorder 
at the begining of the calculation, it cannot give any information about the region of concentrations where the weak 



dilution concept holds 72 



Monte Carlo simulations. Monte Carlo (MC) studies of the d — 3 RIM systems last for almost two decades [|38| pl| . 
One of the first studies of critical behavior of a RIM on a simple cubic lattice revealed the critical exponents in 
a wide dilution region with no deviations within the numerical error from the corresponding exponents of the pure 
(undiluted) system (see Table [il]). However these data were objected by MC simulations on larger lattices p9, ft0[, 
where critical exponents varying continuously with the magnetic sites concentration x were obtained. Indications of 
change of the order parameter critical exponent [3 upon dilution initiated the extension of studies to determine the 
other critical exponents and to check the scaling in the disordered systems. The application of the Swedsen-Wang 
algorithm to the d = 3 RIM |4^] resulted in critical exponents for the susceptibility and correlation length independent 
of concentration over a wide range of dilution. 

Due to Refs. and especially (^H[^6) it became clear that the concentration dependent critical exponents found 

in MC simulations are effective ones, characterizing the approach to the asymptotic region. The effective exponents 
7, j3 and C = 1 — (the last one describes the divergence of the magnetization-energy correlation function) were 
shown p5[ to be concentration dependent in the concentration region 0.5 < x < 1. These data were refined three 
years later f46f| resulting in more accurate estimates for the above mentioned exponents and the critical exponent v 
of the correlation length yielding continuously varying values. The general conclusion of Refs. |^,^] was: while a 
simple crossover between the pure and weakly random fixed point accounts for the behavior of systems above x ~ 0.8, 
in more strongly disordered systems a more refined analysis is needed. Note that this last statement is supported by 



the conjecture of a "step-like" universality of the three dimensional diluted magnets 47 1. 



The critical behavior of the d = 3 RIM was reexamined recently in Refs. |49|-|51|]. In particular, the simulations ]50j ] 
revealed that a disorder realized in a canonical manner (fixing the fraction of magnetic sites) leads to a different 
result from those obtained from a disorder realized in a grand-canonical manner (see Table EB, where the values for 
the second case are denoted by asterisk). The studies of Ref. [^l| were based on the crucial observation that it is 
important to take into account the leading correction-to-scaling term in the infinite volume extrapolation of the MC 
data. The simulations confirmed the universality of the critical exponents of the d — 3 RIM over a wide region of 
concentrations. In particular the value of the correction-to-scaling exponent uj was found to be w = 0.37 ± 0.06 which 
is almost half as large as the corresponding value in the pure d = 3 Ising model to = 0.799 ± 0.011 0. The smallness 
of uj in the dilute case explains its importance for an analysis of the asymptotic critical behavior. 
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III. THE MODEL AND THE RENORMALIZATION GROUP PROCEDURE 



In this section, we define the RG procedure as well as the main quantities which we are going to calculate. Making 
use of the replica method JrjJ and taking the average over different configurations of quenched disorder it is possible 
to show that the RIM critical behavior in the Euclidian space of d = 4 — e dimensions is governed by an effective 
Hamiltonian with two coupling constants [p2| : 

H{ V )= /d^E0 V ^l 2 + *]+| + 1?5>«}' (1) 

J a=l ' \a=l / ' a=l 



in the limit n->0. Here, tp a = ip a (R) is the a's replica of a scalar field; uo ~ x > 0, vq ~ x(x — 1) < is the bare 
coupling constant of the fluctuation's effective interaction due to the presence of impurities; Too is a bare mass. 

In order to describe the long-distance properties of the model (Q) in the vicinity of the phase transition point we shall 
use the field-theoretical RG approach. The results for the RG functions corresponding to ([[]) are obtained on the 
basis of the dimensional regularization and the minimal subtraction scheme |54j|. The renormalized fields, mass and 
couplings <p, to, u, v are introduced by: 
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Here, [i is the external momentum scale and Z$, Z m 2, Za iU , Z^ v are the renormalization constants. They are 
determined by the condition that all poles at e = are removed from the renormalized vertex functions. 

The RG equations are written bearing in mind that the bare vertex functions Tq are calculated with the help of the 
bare Hamiltonian (|l|) as a sum of one-particle irreducible (1PI) diagrams: 

K({r})=<^r 1 )...^(r N )> 1P1 . (2) 



The do not depend on the scale fx, and therefore their derivatives with respect to \i at fixed bare parameters are 
equal to zero. So one gets 

^10 = ^^10 = 0, (3) 



where the index |q means a differentiation at fixed bare parameters. Then the RG equation for the renormalized 



vertex function reads: 



/ d d d d N \ 

K^ +Pu d^ + Pv d~v + lmm dm- " 2 J V ^ K M) = °' (4) 



and the RG functions are given by 



(3 u {u,v) = M^lo) 
<91nZ 

70 = 2-/ 2 {u,v) = n — | , 

ainTO 1 dlnZ^i 

7rn(u,v)=»—\ = -H—^ r \ . 
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Using the method of characteristics the solution of the RG equation may be written formally as: 

r£(m, u, v, M ) = X{£) N ' 2 T%{Y{£)m, u(£),v(£), fi£), (5) 

where the characteristics are the solutions of the ordinary differential equations (flow equations): 

e±\nX(l)= 1 4u(l),v(£)), i^lnY(£) = lm (u(£),v(£)), 

e Ti u{£) = PuW) ' e Te v{i) = PvMQMQ) (6) 

with 

X(1) = Y(1) = 1, u(l) = u, v(l)=v. (7) 

For small values off, the equation (|5|) maps the large length scales (the critical region) to the noncritical point £ = 1. 
In this limit the scale-dependent values of the couplings u(£), v(£) will approach the stable fixed point, provided such 
a fixed point exists. 

The fixed points u* , v* of the differential equations (^) are given by the solutions of the system of equations: 

u (u*, v*) = 0, 

v*) = 0. (8) 

The stable fixed point is defined as the fixed point where the stability matrix 

Bij = , u, = {u, v} (9) 

possess eigenvalues u>x,u>2 with positive real parts. The stable fixed point, which is reached starting from the initial 
values in the limit £ — * 0, corresponds to the critical point of the system. In the limit I — ► (corresponding to the limit 
of an infinite correlation length) the renormalized couplings reach their fixed point values and the critical exponents 
r\ and v of the pair correlation function at T c and of the correlation length respectively are then given by 

i] = 2j 2 (u*,v*), 

l/i/ = 2(l- 7m (iAO). (10) 



In the nonasymptotic region deviations from the power laws with the fixed point values of the critical exponents are 



governed by the correction-to-scaling exponent u> = min(u)x,W2) in accordance with the Wegner expansion 73 



The rest of the critical exponents are obtained by familiar scaling laws. We note, that the expression for correlation 
length exponent may be recast in terms of the renormalization constant Z^i of the two-point vertex function with 4> 2 
insertion by a substitution 2~/ m = 7^+702, which follows from the relations Z m i = Z^iZV 1 and 7^2 = fj,d\nZ~2 /<9/i|o- 



IV. THE RESUMMATION AND THE RESULTS 



In this chapter we analyse series for the RIM f3- functions and critical exponents. The RG functions of the corresponding 
effective Hamiltonian (0) in the replica limit n = have been obtained up to 5-loop order from the appropriate 
expressions for RG functions of the 4 -theory with 0(n)-symmetric and cubic interactions [^3| and read: 
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I3 u /u = -e + 3m + 4m - 17/3m 2 - 46/3mm - 82/9m 2 

+ 32.54968284u 3 + 123.1987313mm 2 + 158.1816418m 2 m + 60.32526811m 3 

- 271.6057842m 4 - 1318. 116311mm 3 - 2452.429994uV - 2003.560971m 3 m 

- 559.7143854m 4 + 2848.568254m 5 + 16789.89843mm 4 + 40367.08593m 2 m 3 

+ 48971. 12730mV + 29091. 77179m 4 m + 6377.751189m 5 , (11) 
(3 v /v = -e + 8/3v + 2m - 14/3m 2 - 22/3mm - 5/3m 2 

+ 25.45714897m 3 + 62.25499170m 2 m + 36.36645522mm 2 + 7m 3 

- 200.9263690m 4 - 667.3761895m 3 m - 650.5641816m 2 m 2 - 259.2586891mm 3 

- 39.91261012m 4 + 2003.976188m 5 + 8469.158907m 4 m + 11721.60876mV 

+ 7434.635066m 2 m 3 + 2344.277996mm 4 + 301.5110976m 5 , (12) 
72 = 1/18m 2 + 1/6mm+ 1/12m 2 - 1/27m 3 - 1/6m 2 m - 3/16mm 2 - 1/16m 3 
+ 125/648m 4 + 125/108m 3 m + 145/72m 2 m 2 + 65/48mm 3 + 65/192m 4 

- 1.005978154m 5 - 7.544836154m 4 m - 18.04854621mV 

- 19.07838990mV - 9.627924878mm 4 - 1.925584976m 5 , (13) 
7 m = 1/3m + 1/2m - 5/18m 2 - 5/6mm - 5/12m 2 + 37/36m 3 + 37/8m 2 m 

+ 251/48mm 2 + 7/4m 3 - 5.3808017m 4 - 32.28481m 3 m - 57.177011mV 

- 39.765731mm 3 - 9.97815253m 4 + 37.850485m 5 

+ 283.878638m 4 m + 686.375317m 3 m 2 + 737.493196mV 

+ 376.1776339mm 4 + 75.37777445m 5 . (14) 



In order to obtain the qualitative characteristics of the RIM critical behavior one can proceed in two ways. The 
historically first scheme is known as e-expansion and consists (i) in expanding the values of the couplings at the 
stable fixed point in e, (ii) inserting these expansions into the field theoretic functions for the exponents (|l3|), ( |l4| ) 
and (iii) expanding again in e. Due to the degeneracy of the /3-functions ( 0), ( |12|) at 1-loop level, an expansion in 
\fe has to be performed |22]-|2^]. Then in 5-loop order on the basis of Eqs. (|ll|) (|14|) one obtains p^,[2^|: 

v = 1/2 + 0.08411582e 1/2 - 0.01663203e + 0.04775351e 3/2 + 0.27258431e 2 , 
r} = -0.00943396e + 0.03494350e 3/2 - 0.04486498e 2 + 0.02157321e 5/2 , (15) 
7 = 1 + 0.16823164e 1/2 - 0.02854708e + 0.07882881e 3/2 + 0.56450490e 2 , 
wi = 2e + 3.704011194e 3/2 + 11.30873837 e 2 , (16) 
u 2 = 0.6729265850 e 1/2 - 1.925509085 e - 0.5725251806 e 3/2 - 13.93125952 e 2 . 



Another method consists in (i) fixing the value of e i. e. the lattice dimensionality, (ii) solving the system of 
equation for the fixed point and (iii) substituting the fixed point values of the couplings into the series for the critical 
exponents |35) (the so-called 3d approach). 

One should note here that often, for the sake of convenience, within the 3d as well as the e expansion one deals 
with the expansions of some combinations of the critical exponents instead of working with them directly on the 
basis of expressions (|l0|). In the present paper the values of critical exponents are calculated from the expansion for 
1/v —1 = 1 — 27,„ and the inverse exponent for magnetic susceptibility 7 _1 = (1 — 7 m )/(l — 72). The numerical 
values of the other exponents are obtained by the familiar scaling laws. 

It is well known that the perturbation theory series for the RG functions in the weak coupling limit as well as in 
e-expansion are asymptotic at best. In order to compare the results obtained on the basis of the *Js expansion and of 
the 3d approach we have to refer to resummation procedures in the calculation of critical exponents. Adjusting the 
resummation procedure we discuss first the one-variable case in both schemes. We start from the £-expansion of the 
pure Ising model critical exponents which in the 5-loop approximation reads JriJ : 

v = 1/2 + l/12e + 0.04320988e 2 - 0.01904337e 3 + 0.07088376e 4 - 0. 21701787s 5 , 

7 = 1 + l/6e + 0.07716049e 2 - 0.04897495e 3 + 0.14357422e 4 - 0.44662483e 5 . (17) 
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An analysis of the e-expansion case starts by representing the expressions for the critical exponents v and 7 of the 
pure Ising model ( p"7| ) in the form of the Pade approximant: [M /N] (x) — Yli=o a i x% I SjLo ^i x ^ m ^ nc variable x = e. 



The results are shown in the form of a Pade table (Table III). The number of the row, N, and of the column, M, 



corresponds to the order of the numerator and the denominator of the Pade approximant [M/N] respectively. One 
can see from this table the expected convergence of the values in the diagonal and first off-diagonal. Therefrom we 
estimate the values of the critical exponents to be v — 0.628,7 = 1.236. These values can be compared with the 



most accurate values v — 0.628 ± 0.001, 7 = 1.234 ± 0.002 |75|, obtained by means of more sophisticated resummation 



procedures |76|] on the basis of the 5-loop e-expansion. We conclude that this good agreement justifies the application 
of the Pade analysis for the e-expansion series (n7h [[77| ■ 



Possessing the information on the asymptotic divergence of the e-expansion we apply a more complicated Pade- 
Borel resummation technique which takes into account the factorial divergence of the series terms. The resummation 
procedure consists in several steps: 

• starting from the initial sum S of L terms one constructs its Borel-Leroy image: 

s = y ad ^y " iixty N , as) 

2—0 I — 

where T(x) is the Euler's gamma function and p is an arbitrary nonnegative number; 

• the Borel-Leroy image ([l8|) is extrapolated by a rational Pade approximant 

[M/N] (xt); 

• the resummed function S res is obtained in the form: 

POO 

S res = dtcxp{-t)t p [M/N](xt). (19) 
Jo 



The values of the critical exponents v and 7 obtained with the Pade-Borel resummation for different M and N are 
given in Table [V. As far as the Pade approximant enters the integral, it might happen that the integrand contains 
poles. If this is the case, the corresponding number in Table IV is calculated by analytic continuation taking the 
principal values of the integral. For the sake of completeness we include such numbers in Table ^ as well as in the 
forthcoming Table VII . However, the final results will be displayed on the basis of data which did not require such 



analytic continuation. Except from the [1/1] case all our final results are obtained from approximants with a linear 
denominator (second columns), which reconstitutes the sign-alternating behavior of the initial series (|l7|). For the 
Borel resummed Pade approximant [4/1] the estimates for the exponents v = 0.629, 7 = 1.236 are in a good agreement 
both with the above Pade analysis (Table III) as well as with the data of Ref. given above. 



In order to complete the study of the pure Ising model we perform an analysis based on the 3d approach. We resum 
the corresponding RG functions of the pure model (they can be obtained by putting v — in the diluted model RG 
functions <JTl|) — <JT^) ) by means of the Pade-Borel resummation technique with a linear denominator approximant (see 
Appendix A) . The results obtained with this method are shown in Table |v|. One should compare them with the 
results obtained recently from the RG-functions in the 3d massive field-theoretical approach, v = 0.6304 ± 0.0013, 
7 = 1.2396 ± 0.0013 [||, and the results of 3d minimal subtraction scheme in the 5-loop approximation v = 0.629 ± 
0.005,7 = 1-235 ±0.005 |55|. Comparing all the results for the pure Ising model we conclude, (i) that the application 
of Pade-Borel resummation technique yields accurate results and, (ii) both the the e-expansion and 3d approach lead 
to reliable values for the exponents. 

We now turn to the results of the Pade and the Pade-Borel resummation techniques applied to the -^/e-expansions ( jlq ) 
( |l6| ) and the 3d approach RG functions (|ll])-([l4|) of the weakly diluted Ising model. We construct Pade approximants 
and perform the Pade-Borel resummation introducing an auxiliary variable t in the expressions (|l^)-(l^) by the 
substitution: e — > et 2 and putting t = 1 in the final results. 



The ap prop riate values for the -\/£-e x P& n sions are listed in Tables VI and VII in the same notations as in the Tables 
III and ttVl One can see from the tables that neither the experimental and nor the reliable theoretical values listed 
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in Section [Il| are obtained. Moreover, considering the expansions for the stability matrix eigenvalues (|16|) it turns out 
that no stable fixed point exists in a strict ^/e-expansion, even with resummation. 

Comparing with the corresponding data for the pure Ising model we conclude that the nature of the y / e-expansion 
series does not allow to obtain reliable information at d = 3 by means of the methods mentioned for the case of pure 
Ising model. This can be considered as an indirect evidence of the nonasymptotic nature of the y^e-expansion. Thus 
a different kind of analysis for the i/e-expansion has to be developed. The fact that e-expansion will not be able to 
give information on critical exponents in system with quenched disorder was predicted already in Refs. pl||6^ ] . There, 
studying the randomly diluted model in zero dimensions, it was shown that the non-Borel summable properties of 
the perturbation theory series are the direct consequence of the existence of Griffiths-like singularities |6S| ] caused by 
the zeroes of the partition function of the pure system. 

In order to treat the theory directly at d = 3 we need a generalisation of the Pade-Borel resummation technique to 
the case of two variables since the RG functions of the weakly diluted Ising model depend on two couplings. The 
corresponding Chisholm-Borel resummation technique can be defined as follows f3~f)] : 



• constructing the Borel-Leroy image of the initial Lth order polynomial S in the variables u and v: 

S= y a, jU V=> V ^i , 

... h. i o<^<L r (*+3+p +i r 



where r(x) is the Euler's gamma function and p is an arbitrary nonnegative number; 

• extrapolating the Borel-Leroy image ( pp| ) by a rational Chisholm approximant [[78j [M/N] (ut, vt) which can be 
defined as a ratio of two polynomials both in variables u and v, of degree M and N so that the first terms of 
its expansion are equal to those of the function which is approximated; 

• the resummed function S res then reads: 

POO 

S res = / Aexp(-i)i p [M/N] (ut,vt). (21) 
Jo 



Here, similarly to the pure case we restrict the approximants to linear denominators and choose the value of the fitting 
parameter p = 0. The motivation of such a choice is discussed in detail in Appendix A. 

Treating the /3-functions by means of this resummation technique leads to a random fixed point, u* > 0, v* < 0, of 
the model already in the 2-loop approximation. The stability analysis shows that this fixed point is stable proving 
the crossover to a new critical regime under dilution. In Figs. [IJ || we show the curves f3 u (u, v) — 0, (3 v (u, v) = in 
the u — v plane. The intersections of these curves (i.e. simultaneous zeros of both /3-functions) correspond to the 
fixed points, the stable and unstable points are marked by open circles and filled boxes respectively. The "naive" 
analysis of the /3-functions, without applying any resummation procedure leads to the curves, which are shown on 
the left-hand side of the Figs. [I], ||. Without resummation only in the 3-loop approximation one gets stable a random 
fixed point u* ^ 0, v* ^ 0. However the fixed point disappears in the 4-loop approximation. A completely different 
picture is observed when the resummation procedure is applied (right hand columns of the figures). In the region 
of interest for the values of the couplings the topology of the fixed point picture remains stable increasing the order 
of approximation from the 2-loop to the 4-loop level; the same behavior has been observed H] for the /3-functions 
obtained in the massive 3d scheme |35| . 

The corresponding values of t he ra ndom fixed point coordinates, critical exponents and eigenvalues of the stability 



matrix are listed in the Table VII] . One can see that increasing the order of approximation one reaches convergent 
results, compatible with experimental and theoretical data (see sec. ||). 



Considering the estimation of the accuracy of the results we note that setting v — in the RG functions (|ll|)-(|14|) they 
are transformed into the appropriate functions of the pure Ising model. In this case the deviation of our 4-loop result 
(see Table 0) from the most accurate one obtained within massive field RG scheme on the basis of 6-loop expansion 
for /3-functions and 7-loop expansion for 7-functions Q is within parts of a percent. This can be considered as a 
upper bound for the accuracy of RIM results. Here the comparison of our data with the results of 4-loop massive 
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scheme results |35],|3(| yields an accuracy of several percents. Since the series for (3- and 7-function (]ll|)-(Jl4|) are 
sign-alternating, the unknown exact stable point coordinates and critical exponents must lie between the 3- and 
4-loop values giving the same estimate (see Appendix A about adjusting a free parameter for the fastest convergence 
of results). Thus, we conclude the accuracy of the RIM critical exponents obtained in our study to be of order of 
several percents. 



A peculiarity of the Table VIIl is that in 5-loop order the applied resummation technique does not lead to a real 
root for the random fixed point. It is expected that the applications of more sophisticated resummation procedures 
incorporating the higher order behavior, still unavailable, will permit an improvement of the estimates for the critical 
exponents in the 2-, 3- and 4-loop level as well as to obtain them in the 5-loop level. However, it is not excluded that 
the absence of a fixed point solution on the 5-loop level might be connected with a (possible) Borel-nonsummability 
of the series under consideration. In this case the 4-loop approximation will be an "optimal truncation" for the 
resummed perturbation theory series, similar to the nonresummed asymptotic series, e.g. in the e-expansion of 0(n)- 
symmetric </> 4 model, where "naive" interpretation of the scries truncated by e 2 term leads to the best (optimal) result. 
On the other hand, numerical and analytic studies of a toy model of a disordered system |3l],[32| revealed two possible 
regimes of the high-order behavior: the first one corresponds to a Borel-summable series whereas the second one does 
not correspond to a Borel-summable series. Numerical studies of up to 200 terms of the expansion resulted in 
the conclusion, that the convergence of the Borel-resummed results depen ds on the strength of disorder (relation of 



the couplings u/v). The convergence of the numerical data of Table VIII is evidence of the fact that the fixed point 
values u*,v* in d = 3 lie in a region, where the resummed series gives reliable information. In any case, on the basis 
of the above analysis one can definitely say: while for the pure Ising model the e-expansion and the 3d approach 
analysis of the RG functions are of equal usefulness, an interpretation of these functions in the diluted model can be 
done only within the framework of the 3d approach. The application of the -y/e-expansion remains to be valid only 
near the upper critical dimension d = 4. This conclusion holds at least within the discussed resummation schemes. 



V. EFFECTIVE CRITICAL BEHAVIOR 



If one is not within the asymptotic region of the stable fixed point power laws for the physical quantities may only 
apply with effective exponents. The critical behavior is then to be understood as a crossover behavior between 
the uncritical background behavior and the true asymptotic behavior. As it has been noted in Ref. |37j| this has in 
principle nothing to do with crossovers between special fixed points (e.g. the pure one and the random one). However, 
depending on the nonuniversal parameters entering the nonasymptotic behavior, the crossover may be more or less 
influenced by the unstable pure fixed point. 

The effective exponents are defined by logarithmic derivatives of corresponding thermodynamical quantities with 
respect to reduced temperature t fl52f . In the RG scheme they are calculated in the region where couplings have not 
yet reached their asymptotic (fixed point) values and depend on the flow parameter I. For instance the magnetic 
susceptibility effective exponent ~f e ff is defined by: 

7e//(r) = = 7(«(*(r)),i;(*(r))) + . . . , (22) 



where the second part is proportional to the /3-functions and comes from the change of the amplitude part of the 
susceptibility. In order to proceed we have to neglect this part since we do not know the amplitude function in the 
same orders as the field-theoretical functions for the exponents. Moreover, the contribution of the amplitude function 
to the crossover seems to be small, at least in other cases [^9|,|8^] . This approximation has also been used in the earlier 
work on diluted models p7j . 

The flow parameter t(r) may be related to the temperature via the matching condition m{l) = (S.q 1 ^) 2 and TOq ~ r. 
However, we discuss the effective exponents as functions of the flow parameter. Then the effective exponents are 
simply given by the expressions for the asymptotic exponents ( |Io| ) but with replacing the fixed point values of the 
couplings u* and v* by the solutions of the flow equations @: 

l/v e ff = 2(1 - J m (u(£),v(£))), Veff = 2 l2 (u(l),v(e)), 7 e// = Veff(2 - V eff) (23) 
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These solutions are shown in Fig. |3| for several initial conditions. There are shown the two unstable fixed points: G 
(Gaussian fixed point) and P (pure Ising fixed point) and the stable random fixed point R, with both couplings non 
zero. In the background region the couplings are expected to be small thus all the flows shown start near the Gaussian 
fixed point with different ratios v(£q)/u(£q) (curves 1 to 6 except curve 2). Curve 1 is the separatrix connecting the 
fixed point G with P, and curve 2 is the separatrix connecting P with R. All curves starting with a negative coupling 
v remain negative but the dependence might be nonmonotonous (see curve 6). Thus several scenarios for the values 
of the effective critical exponents are possible (see Figs. |I] ||). 

Both in experiment as well as in computer simulations (see Tabs. [| ||) exponents reported differ and even exceed 
the known asymptotic theoretical values. This nonuniversal behavior might be related to the possible nonasymptotic 
behavior found in our different flows as has been suggested in |}7|, [l6| . The difference might be due to (i) the different 
temperature regions of the experiment and/or (ii) the different concentrations. The initial values for the couplings 
in the flow equations depend on the value of the concentration, especially for small dilution one expects the coupling 
v to be proportional to the concentration. If this is the case we expect a monotonous increase of the values of the 
effective exponent to the asymptotic value. A typical scenario is seen in curve 3 of Figs. |, @, |. In this case effective 
exponents equal to the pure model critical exponents might be found in relatively wide region of temperature. Then 
as the attraction region of the fixed point P becomes weaker and weaker, an overshooting is possible and effective 
exponents larger its asymptotic values might be found. This scenario is predicted for larger dilutions and represented 
by the curve 6 of Figs, pi 0, 0. Curves 5, 6 correspond to situation when crossover from the mean field behavior 
towards the random one is not influenced by the presence of a pure fixed point. 



VI. CONCLUSIONS 



We studied the critical behavior of the three-dimensional weakly diluted Ising model asymptotically close to the 
critical point, in the intermediate region and far from it. To this end we calculated the values of the asymptotical 
critical exponents, analyzed the behavior of the effective critical exponents and obtained the value of the correction- 
to-scaling exponent entering the Wegner expression which describes the approach of singular quantities to the critical 
temperature. Our study is based on the 5-loop minimal subtraction field-theoretical renormalization group functions 
p3| of the </> 4 -theory with 0(n)-syrnmetric and cubic interactions which in the replica limit n — > correspond to the 
diluted Ising model case. As the minimal subtraction scheme allows to develop either the -^/i-expansion or the 3d 
approach, we adopted the both schemes and compared the results obtained on their basis. 

The perturbation theory RG functions series are asymptotic at best. In order to calculate the critical exponents we 
adopt different resummation procedures within the both schemes, testing them on the well established case of the 
<j) -model with one coupling (pure Ising model). While after the resummation both e-expansion and 3d approach 
give reliable results for the numerical values of the pure model critical exponents, in the case of the diluted Ising 
model the bad convergence properties of the y^^Gxp^si 011 do not allow to obtain reliable values at d — 3. Using 
the direct resummation of the RG functions at d = 3 in the minimal subtraction scheme we obtained numerical 
values for the asymptotic critical exponents of the diluted Ising model and estimated their accuracy to be of order 
of several percents. The results obtained in the 3d approach agree well with other theoretical and recent Monte 
Carlo simulations. Studying the crossover effect we calculated the effective critical exponents and their flows in the 
nonasymptotic region. We observed several scenarios of crossover in the RIM including: monotonous crossover from 
the mean field values of critical exponents to the random ones; existence of a wide temperature region where the 
RIM effective exponents coincide with the asymptotic exponents of the pure Ising model; possible values of effective 
exponents exceeding those of asymptotic ones. 

Though the 3d approach of calculation encountered difficulties in the 5-loop level we guess that the fixed d approach, 
both within the massive |60]] and minimal subtraction |34| , |65| schemes, remains the only reliable way to study critical 
behavior of the model by means of the RG technique. 

Note added in proof. During processing the article several new results appeared in the field. The massive scheme 
4-loop RG functions of the RIM ]35|, [36| were extended to the five-loop pl[ and to the six-loop level |)2| . While 
the traditional analysis of the first via PB resummation allowed to obtain numerical values for the RIM asymptotic 
critical exponents pi]], the method failed for the higher-order functions |Q. However, an application to the functions 
of a refined resummation procedure which treats renormalized couplings of the RG functions asymmetrically 



10 



reconstituted and improved |95[ earliear data for the RIM asymptotic critical exponents. 
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VII. APPENDIX A: THE CHOICE OF THE RESUMMATION PROCEDURE 



In order to obtain reliable quantitative description of the problem under consideration one should adjust precisely 
resummation procedures necessary for analysis of RG functions. Since no information is available about the high order 
behavior of the series for (3- and 7-functions (|Tl|) - ) ( compare with the d = 2 and d = 3 scalar </> 4 -model where the 
Borel summability of the RG functions is proven ]77|, BjJ and the large order behavior is calculated |82|,|8j|), we reject 
all powerful methods implementing such an information (e. g. resummation refined by a conformal mapping, widely 
used in the models of critical phenomena with one coupling [|,84|) and restrict ourselves to the simplest resummation 



techniques which are the generalisation of the Pade-Borel technique (|is|)-(|l9|) to the two variable case. Among them 
we choose that procedure, which, for the different orders of approximation, provides the maximal stability. That 
means, the fastest convergence of the results as well as the maximal similarity of the topological structure of the lines 
defined by the zeros of the /3-functions. 

a) A d = theory 

Let us start tuning the resummation technique by considering the expressions for a toy-model. It is defined by the 
partition function 



Z = ~^/2 / #1 ■ ■ ■ exp(- ^2 $ - u(^2 <t>if - v ^2 1 
n J -°° i=i 1=1 i=i 

which corresponds to the cubic model described by the Hamiltonian (0) in dimension d = 0. One can easily calculate 
the sufficient number of terms representing Z in the form of a series in u and v for arbitrary n. However, the simplest 
case which reproduces the series in two couplings is n = 2 since the case n = is trivial and n — 1 corresponds to a 
series in a single variable u + v. 

We write the series for Z in the two variables u and v as a series in one auxiliary variable r introduced by the 



substitution u — > ut, v — > vt (the so called resolvent series pq]) and apply the Pade-Borel (PB) method (18)-(|19|) 
to the series in r. Here we choose two possibilities of the approximants, one with linear denominator and the other 
one of diagonal type. Moreover we fix the parameter p (|lq) to p = but discuss other choices for the parameter later 
on. The results of this procedure are shown in Figs. f° r ^ ne specific values u = 0.1 and v — —0.01. One notes 

that a stable convergence to the exact value takes place only for resummations with diagonal approximants, while 
the approximants with linear denominator converge only for certain, "optimally truncated" , orders of approximation. 
The larger the value of the variables u and v the less is the order of "optimal truncation" . We believe that processing 
the divergent series (|ll])-(|l4|) in the same manner one may encounter a similar difficulty. 



Analysing the toy-model series by the Chisholm-Borel (CB) scheme (f20|)-(|2lj), one notes that the coefficients in 
the Chisholm approximant are underdetermined. For example, for the L-loop sum an approximant of the linear- 
denominator type [L/l] beeing uniquely defined requires two additional equations. This is the minimal number 
of additional conditions, necessary, to determine any approximant for L = 2 ... 5 except for L = 4, where the 
[3/2] approximant is determined uniquely without additional conditions. Since, generally speaking, a near-diagonal 
Chisholm approximant requires additional conditions the number of which depends on the order of approximation, we 
reject this type of approximants and consider in the following only the [L/l] type approximants. The two additional 
equations have to be symmetric in the variables u and v, otherwise the properties of the symmetry related to these 
variables would depend on the method of calculation. By the substitution v = all the equations describing the 
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critical behavior of the diluted model are converted into the appropriate equations of the pure model. However, if 
a pure model is solved independently, the resummation technique uses the Pade approximant. Thus, the Chisholm 
approximant has to be chosen in such a way that, for either u = or v = one obtains the Pade approximant of 
the one-variable case. This also needs a special choice of the additional conditions. This is achieved by choosing the 
Chisholm approximant [L/l] with the numerator coefficients at u L and v L equal to zero. 

The analysis of the toy-model series by means of this type of CB technique leads to the existence of an "optimally 
truncated" order of approximation similarly to the one shown in Fig. (fn). 



b) A d = 3 theory 



Let us turn back now to the expressions for the RG functions (pT|)-(pT|) at d = 3. Starting from the PB analysis 
one finds that diagonal-type approximants lead to poles on the real axis. This may be due to the fact that such 
approximants do not reconstitute the sign-alternating high-order behavior of the general term of the RG functions, 
which was confirmed in the particular cases n = 2 and n = 3 [ |S6| . Thus for the further discussions we proceed with 
approximants of linear-denominator type. 

Applying the PB method with linear-type approximants we find that the random fixed point for the d = 3 RIM cannot 
be reconstituted already in the 3-loop approximation, while the random fixed point exists even for nonresummed j3- 
functions in this approximation. The picture does not change qualitatively when we try to increase effectively the 
order of polynomial representation for the /3-functions by resumming expressions [1 + r u (u + v)](3 u , [1 + r v (u + v)\(3 v 
with r u , r v taken as fit parameters (compare with [Q). The modified construction of Ref. |37| for the Pade approximant 
with a linear denominator reconstituting the known large order behavior of the one coupling /3-functions, (3 u (u, v = 0) 
and P v (u = 0,v), does not lead to the appearance of the random fixed point either. We conclude, although the PB 
scheme works in 2- and 4-loop approximations, it appears to be unstable in the 3-loop approximation and therefore 
we eliminate it from our consideration. 



The CB method reconstitutes the random fixed point in 2-, 3- and 4-loop approximation, however it fails in the 
5-loop approximation. In order to reestablish the random fixed point we have varied the values of p as well as of 
r u ,r„, but even then we find no region for r u > 0,r v > so that the random fixed point exists in all orders of 
perturbation theory. We also predicted the values of the sixth order contribution for /3-functions (see Refs. [|37|, [88|| ) , 
but the resummation of such pseudo-6-loop functions did not allow to find the random fixed point either. Again a 
modified construction of the Chisholm approximant with the known linear denominator [p7| failed. Comparing the 



behavior of the toy-model series (see Fig. g) and the convergence of the results of the RIM (Tables VIII) we consider 
L = 4 as the "optimal truncation" order within the CB scheme of linear-denominator approximant. 

Once we have chosen the CB method based on Chisholm approximants of [L/l] type as the tool for analysing the 
RIM RG functions we look now for the fastest convergence of the results with increase of the order of approximation 
in number of loops L. To this end we fit the parameters p entering Borcl-Lcroy images (|o|) of RG functions. For the 
RIM we introduce a measure of total deviation between L- and L'-loop results by a function A 2 = (u*' L ~ u*' L ) 2 + 
(v*' L — v*' L ) 2 + (~f*' L — 7* ,L ) 2 + (v*' L — v*' L ) 2 , where the superscripts L(L') indicate a value obtained in L(i')-loop 
approximation. For the pure Ising model an appropriate measure is given by a similar function with v* = 0. We 
adjust now the parameter p to minimize A. The behavior of A(p) is shown in Fig. |^. A minimum of A for the pure 
Ising model is achieved for p = (curves pointed by boxes and triangles in Fig. ||) in both cases L = 4, L' = 3 and 
L = 5, L' = 4. Similar behavior in p is observed for the RIM (curve pointed by crosses) suggesting the choice p = 
as well. 



c) A cubic model 



Since originally the RG functions under consideration (ll)-dlj) were obtained in order to study the critical behavior 
of the cubic model introduced by the effective Hamiltonian (111), we will use it for another test of the resummation 
procedure. In the case of nonzero values of n the critical behavior of this system is governed either by the 0(n)- 
symmetric fixed point for values of n small enough or by the cubic fixed point for n > n c , where n c is the marginal 
value of the order parameter component. The identification of the marginal dimension n c as well as of the critical 
exponents governing the phase transition of the model was recently performed within the e-expansion [p3| . The same 
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model was studied by means of the massive RG approach |3q ] . The numerical value of n c was found to be only slightly 
less than n = 3 leading to values of the critical exponents of the cubic model practically indistinguishable from those of 
the 0(3)-symmetric model. Recent MC simulations strongly suggest n c — 3 although n c < 3 cannot be excluded flS9[ . 
We present estimations of n c obtained from the conditions that the 0(n)-symmetric and the cubic fixed point of 
the resummed /3-functions coincide. This calculation has two advantages: (i) one tests once more the resummation 
methods and, (ii) a new estimate of n c is obtained from the 3d approach within minimal subtraction scheme. Up to 
now n c has been calculated in the frames of e-expansion in the 5-loop approximation |2^,^3| as well as in the massive 
scheme in 4-loop approximation |35| . We perform estimates within the 3d minimal subtraction approach IpEjl in 2-, 

3- and 4-loop approximation. The corresponding values on the basis of /3-functions resummed by the Chisholm-Borel 
method read: n c (2Loops) = 2.7730, n c (3Loops) — 3.1078, nJALoops) — 2.9500. This should be compared with the 
result of the 5-loop e-expansion n c = 2.958 (Pade analysis) J53| |, n c — 2.855 (Pade-Borel resummation) |27| and the 

4- loop massive 3d RG scheme n c — 2.90 (Chisholm-Borel resummation) |3f|. In Fig. ^ we show the lines of zeros for 
the resummed (3— functions as well as the fixed points for different n. One can see that the change of stability of the 
fixed points appears for n very close to n = 3. As our numerical result yields n c < 3 the cubic model at n = 3 is 
governed by a new set of critical exponents which read 7 ~ 1.387, v — 0.709, a — —0.127, r\ = 0.044. The coordinates 
of the stable cubic fixed point u* = 0.0064, v* — 0.3950 and stability matrix eigenvalues io\ = 0.0440, 0J2 — 0.0055 
should be compared with the corresponding values for the nearby situated unstable 0(n)-symmetric fixed point: 
u* = 0, v* = 0.4009, = — 0.0056, W2 = 0.0751. Of course no experiment can distinguish between the critical 
exponents of these almost coinciding fixed points. Note however that at n > n c the stable fixed point can be reached 
only for u > 0. Then n c < 3 means that all systems described by the phenomenological Landau-Ginsburg Hamiltonian 
with cubic anisotropy of negative coupling should undergo a weak first-order phase transition at n = 3 |9C|] . 
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TABLE I. The experimentally measured critical exponents of the materials, which correspond to the d = 3 RIM. The 
measurement procedures as well as materials specifications are given in the following notations: NMR — nuclear magnetic 
resonance; LB — linear birefringence; NS — neutron scattering; MS — Mossbauer spectroscopy; SMXS — synchrotron magnetic 
x-ray scattering; XS — x-ray scattering; DAG — dysprosium aluminium garnet; r min denotes the minimal value of the reduced 
temperature reached in an experiment. 
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0.349±0.008 








HT 


1983 


Fe x Zn 1 _ x F 2 , 
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2 ■ 10 -11 
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1.44±0.06 


0.73±0.03 


-0.09 ±0.03 
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1985 
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+1% Y powder 
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4 • 10~ 2 


0.350±0.01 
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0.73 
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1986 


Fe x Zn 1 - x F 2 , 
x = 0.46 
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1.5 • 10" a 
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1.31±0.03 


0.69±0.01 
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1986 


Fe x Zm- x F 2 , 
x = 0.9925 - 0.95 
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0.36±0.01 
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1986 


Mn x Zni- x F 2 , 
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1988 
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3- 10" 4 
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x = 0.93 
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1.14 ■ 10" 4 




1.34±0.06 
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1G 



TABLE II. The critical exponents of the d — 3 RIM as they are obtained in MC simulations for different values of magnetic 
sites concentration x. Maximal number of lattice sites simulated is N = L 3 . (The asterisk at concentrations denotes that 
disorder was realized in a grand-canonical manner.) 
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17 



V 





1 


2 


3 


4 


5 





0.500 


0.571 


0.609 


0.608 


0.648 


0.514 


i 
i 


A t^nn 

u.oyy 


U.O / o 


A CAQ 


n ci a 
U.D1U 


n con 





o 
2 


0.t>4t> 


U.621 


A CO 1 ) 

0.633 


a cop: 
U.620 








o 
O 


A a*7 
0.09/ 


a co i 
U.631 


0.62 / 


o 








A 

4 


0. /32 


U.62D 














5 


0.431 














o 


7 





1 


2 


3 


4 


5 


n 
u 


l.UUU 


1.10/ 


19/1/1 


i iqc: 


1 


u.oyz 


i 


1.200 


1.310 


1.213 


1.231 
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TABLE III. The Pade table for the values of the correlation length (v) and the susceptibility (7) critical exponents from the 
e-expansion of the pure d = 3 Ising model. Here and in Tables |v], [vi 
to the order of numerator and denominator of a Pade approximant, 
be constructed. 



VII the number of the row and of the column corresponds 
"o" means that corresponding Pade approximant cannot 



V 





1 


2 


3 


4 


5 
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0.584 
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TABLE IV. The results of the Pade-Borel resummation of the correlation leng th ( v ) and the susceptibility (7) critical 



exponents from the e-expansion of the pure d = 
the real part of the corresponding value is given. 



3 Ising model. Here and in Tables 



VII, 



VIII 



the c -superscript denotes that 
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0.031 
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0.131 


0.024 
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0.4958 


1.239 


0.632 
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0.040 


0.756 
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0.4877 


1.246 


0.634 


0.097 


0.036 


0.792 



TABLE V. The results of the application of the Pade-Borel resummation ([L/l]) to the RG-functions of the pure d = 3 Ising 
model. Fixed point coordinate and the critical exponents of the pure d — 3 Ising model obtained by Pade-Borel resummation 
in 3d scheme. 
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1.219 


1.783 


1 
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1.257 


1.101 
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TABLE VI. The Pade table for the values of the correlation length (y) and the susceptibility (7) critical exponent 


s from 


the v^expansi 011 of the d — 3 RIM. 
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0.560 
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U.oa2 
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O 
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0.529 c 
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4 
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1.000 


1.149 
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1.141 c 
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1.141 c 
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1.086 c 
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TABLE VII. The results of the Pade-Borel resummation of the correlation length (y) and the susceptibility (7) critical 
exponents from the y^-expansion of the d = 3 RIM. 



Loop 


u* 
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7 
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V 




UJ2 


2 


0.7886 


-0.1208 


1.308 


0.665 


0.006 


0.032 


0.162 


0.542 


3 


0.6968 


-0.2484 


1.293 


0.654 


0.039 


0.022 


0.430 


0.974 


4 


0.7188 


-0.1697 


1.318 


0.675 


-0.026 


0.049 


0.390 c 


0.390 c 



TABLE VIII. Fixed point coordinates, critical exponents and stability matrix eigenvalues of the d = 3 RIM obtained by 
Chisholm-Borel resummation in 3d scheme. 
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FIG. 1. The lines of zeros of nonresummed (left-hand column) and resummed by the Chisholm-Borel method (right-hand 
column) /3-functions in different orders of the perturbation theory: 1- and 2-loop approximations. Thick line corresponds 
to /3 U = 0, thin line depicts f3 v = 0. One can see the appearance of the random fixed point u* > 0, v* < in the 2-loop 
approximation for the resummed /3-functions. Stable fixed point is shown by an open circle, unstable ones are shown by filled 
boxes. 
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0.8 T v 0.8 v 




FIG. 2. The lines of zeros of nonresummed (left-hand column) and resummed by the Chisholm-Borel method (right-hand 
column) /3-functions in 3- and 4-loop approximations. The notations are the same as in Fig. Close to the random fixed 
point the behavior of the resummed functions remains alike with the increase of the order of approximation. This is not the 
case for nonresummed functions. 
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0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 

FIG. 3. Flow lines for the d = 3 RIM. Fixed points G, P are unstable, fixed point R is the stable one. 
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FIG. 4. Effective exponent v e ff versus logarithm of the flow parameter £ for the flows shown in Fig. 
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FIG. 5. Effective exponent 7 e // versus logarithm of the flow parameter I for the flows shown in Fig. 
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0.94 




0.85 1 1 1 1 1 1 1 

2 4 6 8 lO 12 

Order of approximation 

FIG. 6. The estimate of the toy-model partition function Z at u — 1/10 and v = —1/100 in dependence of the order of 
perturbation theory L in couplings u and v. The application of Pade-Borel resummation with linear denominator approximants 
of type [N/1] (solid line) provides convergence to the exact value (dotted line) only for some first orders of approximation. 
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FIG. 7. The estimate of the toy-model partition function Z at u — 1/10 and v — —1/100 in dependence of the order of 
perturbation theory L in couplings u and v. The application of Pade-Borel resummation with diagonal approximants of type 
[N/N] or [(N + 1)/N] (solid line) provides perfect convergence of the estimate to the exact value (dotted line). 
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o 

O 2 4 6 8 lO 

FIG. 8. The dependence of a measure A of total deviation between 3- and 4-loop pure Ising model results (triangles) as 
well as of 3- and 4-loop RIM results (crosses) on fitting parameter p. One can see the minimum for p = 0. The deviation 
between 5- and 4-loop results (boxes) for pure Ising model confirms p = 0, too. 
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FIG. 9. The lines of zeros of the cubic model /3-functions resummed by the Chisholm-Borel method in 4-loop approximation 
for different n. Thick lines correspond to f3 u = 0, thin lines depict f3 v = 0. The filled boxes and open circles show the positions 
of unstable and stable fixed points respectively. One can see that the crossover to a new behavior appears for values of n very 
close to 3. Our estimate yields n c — 2.950 (see the text). 
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